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Abstract. We consider the effect of modified gravity on the growth of large- 
scale structures at second order in perturbation theory. We show that modified 
gravity models changing the linear growth rate of fluctuations are also bound to 
change, although mildly, the mode coupling amplitude in the density and reduced 
velocity fields. We present explicit formulae which describe this effect. We then 
focus on models of modified gravity involving a scalar field coupled to matter, in 
particular chameleons and dilatons, where it is shown that there exists a transition 
scale around which the existence of an extra scalar degree of freedom induces 
significant changes in the coupling properties of the cosmic fields. We obtain the 
amplitude of this effect for realistic dilaton models at the tree-order level for the 
bispectrum, finding them to be comparable in amplitude to those obtained in the 
DGP model. 



PACS numbers: 



1. Introduction 

The origin of the acceleration of the expansion of the universe [1] [2l [3l [4] is stih 
unknown despite more than ten years of intense experimental and theoretical activities. 
Four broad types of explanations have been put forward so far. The oldest and 
also the boldest one posits that the acceleration of the universe is simply due to a 
cosmological constant [5]. Of course this approach is fraught with difficulties and 
clashes with our present understanding of quantum field theory and the structure of 
quantum corrections. A combination of anthropic arguments and the the existence of 
a stringy landscape ^ have been used to justify the cosmological constant solution to 
the acceleration puzzle. Another possibility relies on the plausible existence of large 
voids in the universe [7]. In a sense, we would have been fooled and the acceleration 
of the universe would simply be the result of local inhomogeneities. This implies 
a strong violation of the Copernican principle which still needs to be tested and 
confirmed (see for instance [8]). A more mundane approach following in the footsteps 
of the early universe infiation era would be to assume that a scalar field leads to the 
recent acceleration era O [TOl [11]. Such models of dark energy run into three major 
difficulties. The first one is a fine tuning issue akin to the cosmological constant 
problem and related to the quantum stability of the scalar potential. Another severe 
problem is the apparent almost coincidence between the beginning of the matter era 
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and the start of the acceleration phase. This requires some type of tuning which may 
be avoided using tracking models of quintessence. Finally, a very serious problem 
arises as the mass of the scalar field now is of the order of the Hubble rate and 
therefore potentially mediates a long range fifth force [12]. Although the smallness of 
the scalar mass is inevitable on large scales, it turns out that locally, where gravity 
has been tested, non-linear effects can prevent the observation of a fifth force. This 
is the case for Galileon models [13] where the Vainshtein mechanism [14] is at play, 
for chameleon models [El [161 113 [H IE] where a thin shell effect is present and 
for dilatons [5D] where a variation of the Damour-Polyakov mechanism [2T] can be 
implemented [22l |23l [24]. Recently a new possibility called the symmetron was 
presented with similar properties to the Damour-Polyakov mechanism: 25]. In these 
cases, the acceleration of the universe and the compatibility with local tests of General 
Relativity is guaranteed. Finally, a fourth type of explanation has been advocated, it 
involves a modification of gravity. Examples of such models are the DGP model [2S] 
coming from extra dimensional physics and the f{R) models [171|551|5ni [3D] ■ It turns 
out that these models of modified gravity are nothing but disguised scalar field models 
with the Vainshtein and the chameleon mechanisms respectively at play. Hence in the 
following we will focus on modifications of gravity involving a scalar field. In particular, 
we will consider the cases of chameleons and dilatons. 

Cosmologically the models we are considering differ very little from a pure A- 
CDM model at least since Big Bang Nucleosynthesis (BBN). The only possibility of 
detecting modified gravity effects in this context is at the perturbative level. At the 
linear level and in a first order approximation, it is known that the main feature of 
these models is the existence of a characteristic scale corresponding to the Compton 
wave length of the scalar field [TT] [3T]. For scales larger than the Compton wave 
length, the growth of structure is identical to the one in A-CDM whereas a regime 
of anomalous growth appears within the Compton scale. We extend these results to 
second order in perturbation theory as a first step towards a more general treatment 
of the renormalisation group type [33l [34] . 

The analysis we carry out applies to general scalar-tensor theories but does not 
cover other modifications of gravity such as the DGP and Galileon models. In the 
DGP case, second order perturbation theory has been studied with modifications of 
General Relativity of similar nature as the ones we will present [35] . In particular, we 
find that the perturbation equations contain a time and scale dependent modification 
of Newton's constant as well as the appearance of extra terms at the second order 
level. A striking difference between DGP and scalar tensor theories is the fact that for 
DGP the scale dependence essentially occurs at large scales whereas for scalar-tensor 
theories this phenomenon is present below the Compton wavelength of the scalar field. 
We will see however that these classes of models bring much more important changes 
to the coupling structure than for instance the clustering quintessence model explored 
in [5S]. 

In the following we first introduce the chameleon and dilaton models model we 
are interested in, their motivations and basic properties, in section [2] In section [3] 
we present the generic properties of second order calculations and how modification 
of gravity can affect the mode coupling amplitudes. The sections [4] and [5] explore in 
more details the consequences of a realistic change of gravity, either due to a change 

X In |32) one can find an earfier attempt to quantify the effect of modified gravity, in the context of 
scalar tensor theories, to higher order couplings. The effects of modified gravity were however limited 
to changes in the background properties. 
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of the linear growth rate or the presence of extra couphng terms. We finally conclude 
and discuss our results in view of future experiments. 



2. Scalar Fields and Modified Gravity 

We are mainly interested here in models of modified gravity involving a scalar field. 
This is the case of the Galileon model, f(R) theories, chameleonic and dilatonic models. 
The last three types will be briefly analysed in this section. 



2.1. Chameleon models 

Dark energy models suffer from the usual smallness of the scalar field mass on 
cosmological scales. Chameleon models overcome this problem by introducing a 
coupling of the scalar field to matter leading to a density dependent mass. In a 
dense environment chameleons have a large mass which leads to the absence of gravity 
violation in the solar system and the laboratory. Cosmologically, the chameleon models 
can be viewed as modified gravity models with a constant coupling to matter /? and 
a scalar field with a time-varying mass m{a). In these models, the bare potential 
appearing in the Lagrangian is modified by the interaction with the environment. 
The action governing the dynamics of the chameleon field (/> is given by 

S= J d^xV^g I - \g^'^d^<t>d,(l^ - V (<^) | 

+ (2.1) 



where g^i, is the metric in the Einstein frame, g^v is the metric in the Jordan frame, 
Mpi = (SttGat)^^/^ = K4 ^ is the reduced Planck mass. It is the Ricci scalar for the 
Einstein-frame metric, and are various matter fields labeled by i. The Einstein- 
frame metric and the Jordan-frame metric are related by the conformal 
rescaling 

~g^, = e^'^^/^^-ffM- ^ A\4>)g^,, (2.2) 

where /3 is a dimensionless constant. In chameleon models this constant is assumed to 
be the same for all types of matter, respecting the weak equivalence principle where 
all material particles feel the same metric. 

The Klein-Gordon equation in the presence of non-relativistic matter defines the 
following effective potential for the scalar field 

V,f,{<j>) = V{<j>)+pe^t'l''-K (2.3) 

where p is the conserved matter density in the Einstein frame. It follows that VeS has 
a minimum at 4>m\n{p) which is solution to 

when the bare potential is of the runaway type as in freezing models of quintessence. 
The mass at the minimum of the effective potential is given by 
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As ^ < and > for successful chameleon models, we find that 

> S/J^i/^ , (2.6) 

where we have assumed that the chameleon energy density is negligible compared to 
p, /90 ^ p. In contrast to usual quintessence models, in most cases the mass of the 
chameleon is much larger than the Hubble rate and up to numerical factors of order 
unity, 

~ 3;3ij2— ^£1^ > iy2 ^^^^ 0<A/pi. (2.7) 

<?>min(/0) 

Thus it is possible to have fluctuations in the chameleon field at scales much smaller 
than the horizon scale. This is of cosmological importance for structure formation. 
When Eq. (12.71) is satisfied, the minimum of the effective potential is an attractor, 
and the chameleon therefore settles down there, and evolves slowly since at least Big 
Bang Nucleosynthesis. In that case, the potential energy of the chameleon is nearly 
constant and the model is equivalent to a A-CDM (Cold Dark Matter) model. 

The dynamics of chameleon models are not equivalent to a A-CDM model at 
the level of perturbations. The first order perturbation theory is well-known with 
an anomalous growth of structures for scales k/a <^ m{a), i.e. within the Compton 
wavelength of the chameleon. We will extend this result to second order shortly. 

Another class of equivalent scalar-tensor theories is given by the so-called f{R) 
models with a Lagrangian 

S^^J d'x^gJiR). (2.8) 

In fact these theories are equivalent to scalar field models with a scalar field (pn 
identified as 

^ = e-2/3«40« (2.9) 
dR 

and /3 = 1/ V6. The potential of these theories is given by 

RfR^f 

VI 

and the coupling function is 

A(0flj = e^''"*". (2.11) 

In this guise, we see that f{R) models must behave like chameleon fields in order to 
evade gravitational constraints |3D]. These models essentially behave like A-CDM in 
the recent past of the universe. Their perturbative properties are crucial to try to 
distinguish them from A-CDM. 

In the following we will focus on chameleon models defined by an inverse power 
law potential 

A 4+n 

n'/')-A4 + ^ + ... (2.12) 

where we neglect higher inverse powers of the chameleon field. The effective potential 
has a minimum at 

(Pmin = I -J-^ I ■ (2.13) 



V{c^B)^ml,^^ (2.10) 
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and the mass of the field at the minimum is given in Eq. (|2.5p . For this model 
(assuming <C Mp\), 

In a cosmological setting, the Hubble parameter H ~ yTp/M-py so that m ^ V^H. 
Hence, analogously to the previous case, the field quickly relaxes to the minimum of 
the potential and sits there for most of the cosmological history, at least since Big 
Bang Nucleosynthesis (BBN). The fact that the scalar field can have a mass much 
larger than H is one of the main differences between chameleon models and usual 
quintessence models. At the minimum of the potential, 

K.ff(0mi„)« A^+p. (2.15) 

This implies that the model behaves like a pure cosmological constant in the recent 
past. The matter contribution is simply the usual energy density entering the 
Friedmann equation. 

Using the fact that these models behave like A-CDM, we have during the matter 
epoch p = Pcq(acq/a)^- This implies that the field and the mass have a time 
dependence given by, 

0min ~ a3/("+i) (2.16) 

and 

m2(a) ^a-3("+2)/("+i). (2.17) 
It is important to notice that the comoving Compton wavelength scales like 

m(a)a ^ a-(«+4)/2(n+i) (2.18) 

which decreases with time. Hence a given scale k enters the Compton wavelength at 
late time implying that gravity is essentially modified at short distance in the recent 
past of the universe. 

2.2. General scalar-tensor theories and dilaton models 

So far we have focused on models where the mass of the scalar field is time-dependent 
and the coupling /3 is constant. More general models where the coupling is a field 
function itself can be easily constructed. In the following we consider models with the 
following action 

d^x^gC^{ij'i,A\<i>)g^,), (2.19) 

Notice that fc^((/)) = l/2Mpj leads to a canonically normalised field. A particularly 
interesting case concerns the runaway dilaton of string cosmology in the strong 
coupling regime of string theory [5D]. In this dilatonic case we have 

V{(j)) = A^cj^W^e-^ (2.20) 

in the large (j) limit corresponding to the large string coupling limit and 

k\^)^W\4>)-A\^)^^ (2.21) 

ic^ 
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where 



dhiA 



(2.22) 



(2.23) 



is the field dependent couphng constant, and we have 

In the strong coupling limit of string theory, Damour and Polyakov in |21) have argued 
that the coupling should vanish at a finite/infinite value of 0: 



(2.24) 

which has a minimum at (f>o. The coefhcient A2 ^ Mp^/M^ where is the string 
scale must be large enough to satisfy the solar system tests [23]. 
In these models the effective potential is 

V,sW^V{cj))+A{(b)pm (2.25) 

where pm is the conserved matter density in the Einstein frame. We define the 
normalised field 

dip = k{<j))d<j) (2.26) 



whose mass is 



dip^ 



(2.27) 



where (j) is the cosmological background value. At the background level, like chameleon 
models, the dilaton models behave like a A-CDM model. At the perturbative level, 
we can define a critical wave mode 

k, = air,o)m{4>{r]o)) (2.28) 
below which gravity is modified now. 




Figure 1. The effective mass of the field </3 as a function of time rj in units of 
Hq. The numerical application corresponds to the model H2.24|l with A = 1 (left 
panel), A = 10 (middle panel) and A = 1000 (right panel), = 5.6 10^ (solid line) 
and the dashed lines to cases where A2 is 1/3 or 3 times larger. In this case the 
value of kc is kc = 4:10^ Hq/c 1.3hMpc~^ . Gravity is modified for modes that 
are above the solid line at time rj. 



3. Parametrised Modified Gravity 

In this section we consider a simple parametrisation of modified gravity and explore 
the consequences of this very phenomenological approach. It has been argued in the 
literature that models of modified gravity could be characterized by the growth rate 
of the perturbations. We envisage here the phenomenological consequences of such an 
effect on the higher order correlation functions assuming the change of growth rate is 
entirely due to a (time dependent) modification of the gravitational force. 
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3.1. The ^- model 

As a first approximation before considering more general theories, we focus on models 
where matter is conserved: 

^^ + i[(l + ^)u,(x,i)],i = 0. (3.1) 

and (5(x, t) is the density contrast at comoving position x and physical time t. This 
equation is valid on sub-horizon scales. Here and in the following wc do not make the 
distinction between the dark matter fluid and the baryon fluid. Then in the single 
flow approximation, neglecting pressure terms, the Euler equation takes the form, 

^^^^ + i/(t)u,(x,0 + ^u,(x,t)u,,,(x,i) = --L$f (x,t) (3.2) 

where <E''^'^' (x, f) is the effective potential seen by the massive particles. Before 
exploring more realistic settings, we simply assume here that it is given by 

$°ff(x,t) = (l + e(i))$(x,t) (3.3) 

where $(x, t) is the standard Newtonian potential; i.e. it satisfles the standard Poisson 
equation, 

^ A$(x, t) = AnGpit) 5(x, t). (3.4) 

One can then introduce 9 = Ui^i/{aH) which is nothing but the local expansion rate 
in units of H. Then the conservation equations can be written in Fourier space, 

5{k, t) + He{k, t) = -if 5d (k - ki - ks) a(ki , ka)^ (fci , t) ^(fcs , t) (3.5) 

e{k, t) + i^H + I j e{k, t) + ^H^Qmil + e{t))S{k, t) = 

-ff(5D(k-ki -k2)/3(kl,k2)^^(fcl,^)^^(fc2,^) (3.6) 

where a(ki,k2) = k2.(ki + k2)/A;| and /3(ki,k2) = (ki.k2)(k2.(ki + k2))/fc^/A;| and 
where a dot stands for a derivative with respect to time. The right hand side terms are 
implicitly assumed to be integrated over ki and k2. We have furthermore introduced 
the density parameter fim(i) which represents the (pressureless) matter density in 
units of the critical density. 

When linearized this system leads to a scale independent linear system for the 
growth rate D, 

bit) + 2Hb(t) - ^H'^Q.^{1 + e{t))D{t) = (3.7) 

Let us denote by _D_|_ the solution of this equation which grows with time (assuming 
there is only one of such solutions). This is the case for standard gravity, for which 
e = 0, so by continuity it should be similar when e is small enough. We can then 
define the factor f{t) as 

_ d\ogD+{k,t) 

A\oga{t) ■ ^^-^^ 

so that the time dependence of 9{k, t) for the growing mode is nothing but f{t) D+{t). 
In this section, and in this section only, we assume that f{t) takes the form, 

m = ^i, (3.9) 
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and that e(i) is such that f(t) has such an expression when rini(0 given by the 
background evolution of the standard A— CDM model. It is known (and we will 
recover this result) that for standard gravity 

7°^ w 0.55. (3.10) 

For such a standard background evolution we have H/H^ = — 3/2ilni and Ctm = 
3i?Slni(f^m — 1) so that one can encapsulate the entire time dependence in the flm{t) 
dependence. 

Assuming Eq. (j3.9p . the linearized conservation equations can be used to infer 
the form of e{t): 

1 + ^-rit) = [^'J + sjnun,,, - 1) + (2 - 3n,j2)ni] . (s.ii) 

In the next paragraph we will analyse the consequences of such a change for the growth 
of higher order correlation functions. 



3.2. Bispectra in the "f— model 

The bispectra will be driven, at large enough scales, by the mode coupling evolution 
deduced from the equations of motion. Let us assume that one can expand the cosmic 
fields with respect to the initial (linear) fields, 

5(k) = ^<5(^)(k), e{\.) - ^0(^)(k). (3.12) 

p p 

where (5'^^^(k) is of order p in (5'^^(k). Generically we can write, 

6'^P\\<i) ^ Fp\5^^\\<i)Y (3.13) 

and similarly for the velocity divergence with functions Gp (we will make these forms 
more precise in the following). Then for Gaussian initial conditions the tree-order 
density bispectra are directly related to the second order density field. More precisely 
if we define the density spectrum P{k) and bispectrum B{ki, k2, k^) as 

(5(ki)5(k2)) ={27rfSD{k,+-k2)P{ki) (3.14) 

(5(ki)5(k2)<5(k3)) = {27rfSD{k, + k2 + k3)Bs{ki, k2, fca) (3.15) 
then at tree order 

Bsikuk2,k3) = 2F2{kuk2)P{ki)Pik2) + sym. (3.16) 
Note that for convenience one usually defines the reduced bispectrum as 

D fi 1 I, \ Bs{ki,k2,k3) 

"^'^''^ = P[k,)P[k2)+P[k2)P{k,)+Pik,)P[k,) ^'-''^ 

which is then independent of the amplitude of the spectrum. 

The quantities that are accessible to observations are generically the density field, 
for instance the projected density field in case of the cosmic shear, or, for redshift space 
observations, a combination of the density field and the reduced velocity divergence, 
where the reduced velocity divergence 9{x, t) is defined by 6{x, t) = 0(x, t)/ f{t). In the 
following we simply assume that the bispectrum Qs or similarly that of the reduced 
divergence Qe are both accessible to observations and we will therefore analyse the 
impact of modified gravity models on these quantities. 
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In Fourier space the equations of motion then become, 

5'{k) + 9{k) = -5u(k - ki - k2)a(ki, ka) S{ki) ^(fca) (3.18) 



-(5D(k - ki - k2)/3(ki, k2)^(fci)0~(fc2) (3.19) 

where a prime denotes a derivation with respect to the logarithm of the growing mode, 

FIX 1 

X' = D+^- = ——X. (3.20) 

+ dD+ f{t)H ^ ' 

The system p.l8ll3.19|) can be solved recursively when one assumes that the fields 5 
and can be expanded with respect to the linear solution, 

5(k) = 5]-5(^')(k), 0» = ^0~(^)(k). (3.21) 

p p 

The leading time dependence of S'^'p'^ (and 9'^^'^) can easily be deduced and goes like 
-D^. On can then formally define the functions Fp and Gp as 

6^P\\<i)^ J d3ki...d%fci,ac(k-;^k0 

i 

X Fp(ki, . . . ,kp)5(i)(ki) . . .(5(i)(kp) (3.22) 

^(^^ (k) = - /" d^ki . . . d^kj, 5Dirac(k - J2 

X Gp(ki,...,kp),5(i)(ki)...(5(i)(kp). (3.23) 
In general this system can be solved when ^(i) defined by, 

a, 
p 

is independent of time. The functions Fp and Gp are then independent of time and 
the system p. 18113.1^ gives recursive relations (this extends the relations (43-44) of 
the review paper [57] obtained for an Einstein-de Sitter background in GR) for the 
functions Fp and Gp, 

pi^p(ki,...,kp)-Gp(ki,...,kp) = 

p-i 

^a(qi,q2)Fp_g(kg+i,. . . ,kp)Gg(ki, . . . ,kq) (3.25) 

q=l 

l-p-y)Gp(ki,...,kp) + ^Fp(ki,...,kp) = 

p-i 

- E /3(qi, q2)Gp-,(k,+i, . . . , kp)G,(ki, . . . , k,) (3.26) 

9=1 

with Fi = Gi = 1, qi — 12 = J2^=q+i ^i- When the coefficient is time 

dependent the situation is slightly more complicated but the general forms of the 
functions F2 and G2 is still very constrained. From the continuity equation it is easy 
to see that we have the relation, 

G2(k„k2) = 2F2(k„k2) - -(^^^^^)+ -(^^^^^) (3.27) 



e(i) = ^(l + 6), (3.24) 
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It is slightly more cumbersome to see that at late time the functions F2 and G2 
necessarily take the forn|§|, 



i^2(ki,k2)= -r-- 



T 



and 



3/i2 

4 



2 kj 
1 ki.ka 



2 A;| 

1 ki.ks 

2 fc| 



4 



(ki.k2)2 
(ki.k2)^ 



(3.28) 



(3.29) 



where 1^2 and /i2 are the only time dependent coefficients [Jj which depend on the 
background through Oin(l + e)//^- For Einstein-de Sitter in GR, this is simply unity 
and it is easy to find that 1/2 = 34/21 and /i2 = 26/21. The expressions of 1/2 and /12 
can easily be computed numerically for the 7— model. The results are presented on 
Figs. H 




0.5 0.55 0.6 0.65 0.7 0.75 0.8 

7 




0.5 0.55 0.6 0.65 0.7 0.75 O.S 

7 



Figure 2. Value of U2 (left panel) and fi2 (right panel) as a function of 7 in 
the 7— model for z = (solid line) and z = 1 (dashed line) and the expressions 
l|3. 3613. 3711 (corresponding dotted lines). 



In the standard General Relativity (GR) regime and for a matter dominated 
universe, ^ ~ 1 so that the results must be close to the Einstein-de Sitter case at early 
time (high redshift) and gradually depart from it at late time. To a large extent, the 
results obtained here can actually be derived from a perturbation calculation when 
« 1. Let us assume that 

770 = 1 - (3.30) 

is a small (time) dependent variable. Then from Eq. p. lip one has. 



^7 ~ 2 ( 1 - y7 ) 



(3.31) 



The parameters 1^2 and fi2 can similarly be computed perturbatively about the 
Einstein-de Sitter solution, 



and one finds 



34 ^ 26 , 

^2 = [1 + Vi') , M2 = ^(1 + ??m) 



1 5 



221 



(3.32) 
(3.33) 



§ This comes from the decomposition of the right hand side of Eqs. into the growing and decaying 
modes of the linear system. It turns out that ki.k2/fcj + ki. 142/^2 is always in the growing mode. 
II Actually 1/2 and fi2 can be interpreted as the geometrical average of F2 and G2 respectively and 
more generally the background dependence of Fp and Gp are entirely encoded into their geometrical 
averages which in turns can be derived from the spherical collapse dynamics. See section 2.4 of I37| 
for details. 
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Note that in this approach the GR resuh corresponds to 7 = 6/11 (indeed close to 
0.55) for which vanishes. In this case we actually find that the 1/2 and /i2 parameters 
are well parametrisec0 by 

-r^l + ^^n^^'''- (3.34) 

M^"- -| + y^n^/^^^ (3.35) 
More generally we find that the relations 

- 7^^)(l - ^MT-' (3.36) 

= 7°^)(1 - ^^m)^^^-^ (3.37) 

give a good fit to the numerical results for realistic values of fim. 
4. Linear Scalar- Tensor theories 

In this section we will deal with models where the coupling j3 is constant and the 
potential is the one of a massive particle with a time dependent mass. This is a first 
order approximation to the behaviour of chameleon and dilaton models. The effects 
of the perturbations of the mass term and the coupling function will be dealt with 
later. 

We need to amend the perturbation equations to take into account the effects 
of modified gravity. We first consider models where matter is conserved and matter 
particles interact with a new scalar force of constant coupling and a time varying 
range. The continuity equation is therefore left unchanged. This simply expresses 
that matter (here either baryonic or dark matter) is assumed to be conserved during 
the cosmic evolution. In the single flow approximation, neglecting pressure terms, but 
in the presence of an extra scalar field, the Euler equation takes the form, 

^"'g^'^^ + H{t)n,{^, t) + (x, t)u,j (x, = ('i'(x, t)+p 0(x, t)) ^ (4.1) 

where $ is the usual Newton potential (the one that enters in the metric) and (/)(x, t) 
an extra scalar field that couples to the matter field. Both fields are coupled to the 
matter fluctuations. The potential $ follows the usual Poisson equation Eq. (13. 4p . 
whereas the held obeys the following equation 

^ A0(x, t) - m^c^ = 87rG/3p(t) <5(x, t). (4.2) 
a (i) 

where /3 is a free coupling and A is the flat space Laplacian. It determines the strength 
of the modification of gravity. The other relevant parameter is the mass scale m. This 
mass scale is a priori time dependent We will see later that the fact that ma is 
decreasing with time implies that, for any scale, gravity was always standard in the 
past. 

^ This form for U2 is compatible with the expansion we found and suggested by the expression of the 
second order displacement field as shown in |37| . The expression for ^2 in this reference is however 
incorrect. 

+ In the numerical applications we do here we assume that ma ~ a^'^/^. which corresponds to 
n = 1/2 in terms of chameleon models. 
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This system of equations form a cfosed system between the focal density contrast, 
the pecuhar vefocity u, the gravitational potential $ and (/). Note that these equations 
are valid in the sub-horizon limit, at linear order in (j) (higher order corrections will 
be dealt in the following section on general scalar-tensor theories) and any order in 
the density contrast. As for the standard equations of motions, we can note that the 
source term of the Euler equation is a potential term. As a result one expects the 
peculiar velocity field to derive from a potential too if it was so initially. In Fourier 
space it is then possible to get the very same system as Eqs. ()3.5ll3.6p but changing 
e{t) to a k dependent function that reads 

1 + rn?a^ jk'^ 

which captures the whole effects of modified gravity. As a consequence the linear 
growth rate is now k dependent. 



<fc^^)= . , J„2/,2 (4-3) 



Jf..l. The linear solution 



2. 




0.2 0.4 0.6 0.8 1. 1.2 



Figure 3. Dependence of the function /fe, defined in Eq. 114.51 1 as a function of 
Qm as given by Eq. ||4.10| I. The solid grey lines show the behaviour of for 
standard gravity (lower curve) and modified gravity (upper curve) for 2/3^ = 1 
assuming the the wave mode k is such that k ^ mao. The black lines show the 
actual behaviour of fk{^) as a function of Vlm{ci) following its variation during 
the expansion of the universe for k/{mao) = 100, 10, 1 and 0.1 from the upper 
curve to the lower curve. Depending on k, the transition from the /3 = to the 
13 = l/\/2 solution takes place at different time. 



When linearised this system leads to the following equation for the growth rate 
of the fluctuations, 

D{k, t) + 2Hb{k, t) - ^H^n^^{l + e(fc, t))D{k, t) = . (4.4) 

Let us denote by the solution of this equation which grows with time. We can 
then define rjk = log D+{k,t) and the factor fk{t) as 

d\ogD+ik,t) 

h = — T-j 77^ — ■ (4.5) 

dloga(i) 

Thus we have, for any time dependent function X, 

X = fkHX', (4.6) 
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where X' is the derivative of X with respect to ?7fc. We note that the Unearized 
continuity equation reads 9{k) ~ —fkS{k) and the equation (j4.4l) can be rewritten as, 

2 + h-l^{l + e{k)) + fl, + ^ = Q. (4.7) 

The resulting behaviour of fk is shown in Fig. [3] assuming H / is the one of a 
A— CDM universe, i.e. H/H^ = — 3/2ilni. In this case one has, 

fkn[^ = 3r!m(rim - 1). (4.8) 

The evolution equations for fk and i7m then lead to, 

dfk ^ 3/2an(/ + l + £)-/fc(2 + A.) 

with the initial condition that fk is constant when — > 1 . This equation can actually 
be solved explicitly for fixed values of e. The solution reads. 



42^1 ( /±, i (/+ + 2) ; i (4/+ + 7) ; 1 - ji- ) (4/+ + 7) f}„ 



12 2^^1 



I if + + 3), I (/+ + 5) ; i (4/+ + 13) ; 1 - ^ 



(e + /+ + 1) (^n - 1) 



with 



42F, I y , J (/+ + 2) i i (4/+ + 7) ; 1 - ^ I (6f + 5/+ + 6)n„} (410) 

corresponding to the Einstein-de Sitter value of / in the presence of modified gravity. 

This solution is illustrated in Fig. [3l The two grey lines correspond to constant 
values of e, either e = (which corresponds to standard GR, lower line) or e = 1 
(small scale modified gravity with 2/3^ = 1, upper line). It is found that the general 
solution, for a constant value of e is very close to[f], 

2(2 + /+) 

f{nm)^f+nj^. (4.12) 

When k is comparable to or larger than kc = m oq (ao is the current value of the 
expansion factor) the function fk switches for one curve to the other as shown in Fig 
[3] with the solid black curves. It is to be noted that this result does not follow the 
form of Eq. 13.91 a feature already noticed in [21] . 

4-2. The mode coupling evolution 

Similarly to the first section, let us define the reduced velocity mode by dividing it 
by fkit) thus introducing 9(k) = 9{k)/ fk- The continuity and Euler equations now 
become 

6'{k) + 9{k) = -,5D(k - ki - k2)a(ki,k2) (5(fci) 9{k2)fy- (4.13) 

fk 

^'(k) - - + m)) 0{k) + ?%(1 + eik))S{k) = 



2 Pk ' ' "J ' ' 2 fl 



. fki fk2 



~S^{k - ki - k2)/3(ki, k2)0(fci)0(fc2)^^ip{4.14) 

Jk 

* This form can be obtained assuming /+ f!" and expanding Eq. ||4.9|I to first order in uj with 

n = 1 + a; which gives the value of u. 
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The evolution of the nonhnear couphngs wiU then depend on + e(/c)). 

In general this system can be solved explicitly in a regime where all the /^'s are 
independent on k and where i^m/ fk (1 + ^(^)) independent on both k and time. 

In the standard GR regime and when the universe is still dominated by its matter 
content we have flm/ ~ 1- This is the regime we have at early time, that is at large 
enough redshift. In the context of standard cosmology, it happens that, for a A-CDM 
background, the solution stays close to / = 17"^^ so that F2 and all higher order 
functions Fp are only weakly evolving with time. They remain close to the Einstein- 
de Sitter solution. In the context of modified gravity regime this will still be the case 
for large enough scales (i.e. k ^ kc) since then e{k,t) is small (and has always been 
so assuming that ma is decreasing with time). As a result the higher order terms, 
to all orders, are expected to remain formally the same (as for the Einstein-de Sitter 
case) when they are expressed in terms of the linear solution. This will be illustrated 
on the left hand plateau of Fig. HI 

Another regime of interest is when all scales are well within the modified gravity 
regime. In this case e(fc, t) is both independent of k and of time (it is simply 2/3) and, 
according to Eq. (|4.12|) . rim//| (1 + e(fc)) is also found to be roughly constant but 
with a slightly different value. As a consequence the formal expression of F2 show in 
Eq. p.28p is still valid provided V2 is adequately calculated. In the fim — > 1 limit we 
have 

'4±^ withC=^, 
3(4 + 30 ^ fl 

which incidentally varies from 34/21 to 14/9 when /3 varies from to 00. Quite 
surprisingly 1^2 remains finite even when the effect of self-gravity is made arbitrarily 
large. This induces a change, although modest, for the tree-order bispectrum as we 
will see later. 



^2(e) = :::To.' ^1*^ ^ = (4.15) 




-2-10 1 2 3 4 

Logjo(k/k<.) 



Figure 4. Amplitude of the reduced bispectrum of the density field for equilateral 
configurations as a function of k/kc for an Einstein-de Sitter background (dashed 
line) and a A— CDM background (solid line). For low values of k, the result 
expected for standard GR is recovered. For high values of k, modified gravity 
effects take place. 



In general however one is interested in mode coupling effects for modes that may 
not be in the same regimes and for which we do not have necessarily fk = /^^ = /^^ . 
The general form F2 is then changed and a case by case numerical integration is then 
necessary (see [55]). 
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4..3. Bispectra 

The effect of modified gravity is illustrated in Fig. 2] which gives the expression of the 
reduced bispectrum in case of equilateral configurations (e.g. ki — k2 = fcs). Note 
that in this case, the result is independent of the shape of the power spectrum. We can 
see the effect of modified gravity compared to the effects of the background evolution 
(dashed to sohd line). For the dashed line the asymptotic plateaus can be obtained 
from the form (|3.28|) and the values of V2 given by Eq. (|4.15|) for e = and e = 1. 




Figure 5. Amplitude of the reduced bispectrum Q{ki, k2, fca) of the density field 
for the modified gravity divided by the expected result for standard gravity as 
a function of k\/k2 and k^/k2 for k2 = lOfcc- The left panel corresponds to 
P{k) ~ k~^ and the right panel to P{k) ~ k~^ . Note that the irregular patterns 
on the right sides of the figures are due to inadequate interpolation effects. 

Finally in Fig. [5] a more general k dependence of the bispectrum is shown and 
compared to the standard GR results. It shows that for some elongated configurations 
modified gravity can depart from standard GR by as much as 15%. Notice that the 
effect is expected to be much larger for the statistical properties of 9 as we will see 
it in a more complete analysis. The computations have been done for two different 
hypotheses for the power spectrum shape, P{k) ^ k^^ (left panel) and P{k) ^ fc^^ 
(right panel). The two panels exhibit subtle differences that are due the fact that 
the different contribution to the bispectrum are weighted slightly differently for non- 
equilateral configurations. 

5. Modified Gravity with a General Scalar- Tensor Theory 

5.L Modified gravity 

For a general scalar-tensor theory with a field dependent coupling (3, the perturbation 
equations are slightly modified from the ones used in the linear case. First of all the 
mass is sensitive to the perturbations and therefore the third derivative of the potential 
plays a role. Similarly, the coupling can be perturbed and thence its derivative is also 
relevant. First of all in the Einstein frame the Einstein equations read 

Ft 

G^i, = R^i, — —gfj,v ~ i^i{T^^ + T'^u) (5.1) 
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where T™ is the matter energy momentum and T*'^ is the scalar field energy 
momentum tensor 

T;f, = 2d^^d,.^ - (V{^) + {d^f)g^u (5.2) 

where we are working with the normalised field f for simplicity. The Klein-Gordon 
equation is simply 

V'^V,,^ = ^{d^V - p,s{v)T^) (5.3) 

where T™ is the trace of the energy momentum tensor and /3eff = l3{'P)/k{4>) Using 
the Bianchi identity V^G'^'^ = we find the non-conservation equation 

V^Tr=/3cff(^)T™9V. (5.4) 

In the case of a pressure-less fluid with 

Tr - p'^u'^u-^ (5.5) 

where = dx'^/dr and r is the proper time along the evolution of CDM particles 
with = — 1, we find that 

+ = p^/3eff(v)0 (5.6) 
where F = u^^Vf^F for any tensor and the local Hubble rate is 3h = V^u^. Similarly 
the Euler equation becomes 

+ Piip)ipu^ = -Pes{^)d''ip . (5.7) 

In the following we will neglect the time variation of the scalar field compared to the 
spatial gradients. Hence in the non-relativistic limit the Euler equation becomes, 

e = -/3eff((^)aV, (5.8) 
implying that the CDM particles feel the presence of a scalar force of coupling strength 
l^esif)- It will be very useful to redefine the matter density 

p^ = A{v)p (5.9) 

so that the conservation of matter equation is satisfied (notice that a dot here stands 
for the derivative with respect to proper time) 

p + 3hp = (5.10) 

or equivalently 

V^(pw'') = 0. (5.11) 

This is the usual conservation of matter in cosmology. The conserved matter density 
is not p^ but a rescaled density depending on the Weyl factor A{(fi). 
Finally the Klein-Gordon equation can be written as 

VV^V' = ^ttGnO^VM, (5.12) 
where the effective potential is simply 

VM^) = V{cl>)+A{(b)p (5.13) 

and the conserved density p is ^-independent. This is the expression we have used 
already for chameleon and dilaton models. These equations lead to the proper set 
of perturbation equations in the general case when complemented with the Poisson 
equation with the Newtonian potential depending on A{tp)p. 

In conclusion, we have found that the Euler equation is modified by the presence 
of a scalar force while conservation of matter is still valid. The conserved matter 
density must be dressed with the conformal factor A{ip) to act as a source for the 
Newton potential in the Poisson equation. 
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5.2. Perturbation equations 

The Euler equation can be written explicitly in the single flow approximation, 
neglecting pressure terms: 

-^^-L^ + H{t)ui{x,t) + — u,(x,i)uij(x,t) = 

(a>(x, t),, + /3eff(^ + S^M^, t),,) (5.14) 

where $ is the usual Newton potential (the one that enters in the metric) and (f>{x, t) 
the scalar field that couples to the matter field. Both fields are coupled to the matter 
fluctuations. For the Newton potential we have the usual Poisson equation, 

^ A$(x, t) = 4nGA{4>)p{t) 5m (x, t) (5.15) 

a (t) 

where 

5m(x,t) =<5(A(<^)p) (5.16) 

which depends on both the scalar and matter perturbations. Using the Klein-Gordon 
equation 

1 



-A<p(x,t) + 47rG [V,s,^{^) - V,s,^i4>)] 



a^ity 

AnG p,s{<P)A{4>)m <5m(x, t). (5.17) 

where we have defined 

M) = f| (5.18) 

we find that the Klein-Gordon equation, the Poisson equation and the Euler equation 
involve Sm and not 5p/p, i.e. the Einstein frame energy density. As a result the 
perturbation variables can be conveniently chosen to be Sip, 5^ and 9. Of course this 
is not true of the matter conservation equation as the Einstein frame energy density 
is not conserved in general. On the other hand here the Klein-Gordon equation leads 
to the estimate for the p—th order in perturbation theory 

(5.19) 

implying that the discrepancy between 6p/p and Sm is negligible when H. This 
condition is satisfled for the models we consider as the background mimics A-CDM 
and the background value of the scalar field is an attractor. In the following, we will 
use the perturbed equations as a function of 

5.3. Second order perturbations 

We now focus our attention on the second order perturbation equations. The first 
step is to take the the divergence of the Euler equation: 

{aiie) + aH^e + jV,{ujV,Uj) = -^(A$ + /3eff(0)A5^ + Vi0^s{<t>Wi5^). (5.20) 
Working to second order only, this equation simplifies and becomes 

^AaHe) + aH^e + -ViiujWiUj) + -(A$ + /3eff(^)A<5^) 

do Qi CI 

= -^M4>)^i{SipViSip)) (5.21) 
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^''^^^^ = ^gr- ^^-''^ 

Let us now analyse the perturbed Klein-Gordon equation at the second order: 

A - 3 

— + m^{(j)))Sip+-i 



-^ir2l7^A(^)7eff (0)Mm - (5.23) 

where 

u(0)=47rG4^. (5.24) 

We then get the Fourier components of the first and second order perturbations of the 
field. At first order we have, 

5^(1) = '^S{k)5^^ (5.25) 
with the function S{k) defined as 

The behaviour of this function is depicted in Fig. [SI At second order the Klein-Gordon 
equation reads 

V^(^)=i?^<5(?^-^l]™i?2A(0)7cff(^)- 



(5.27) 



2 H-a2m2((/.)/fc2 
where we have used the notation for the product of perturbations 

AB = fc(k-ki -k2)^(fci)B(fc2). (5.28) 




Log,„(k/kc) Lo8,„(k/ke) Log,„(k/k,) 

Figure 6. Dependence on k of the parameters S{k), r]{k) and fJ-{k) for r^fj = 
(solid lines), rja = —1 (long dashed) and r;n = —2 (short dashed). 

We are now in position to write the Euler equation at second order. It takes the 
form: 

1^(2) + (2 + -^)0(^) + + e{k))S(^ = - /?(ki , k2)0^ 

- [5Eul.(kl,k2) +5lntr.(kl,k2)] {S^^^f (5.29) 
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Li^iiiWIic) LoB,„(k/k.) L()g,„(k/k^) 

Figure 7. The density (top row) and reduced velocity divergence (bottom row) 
bispectra for equilateral configurations as a function of scales for A = 1 left panels, 
A = 10, middle panels and A = 1000, right panels. The solid line is the General 
Relativity prediction. The modified gravity model used here corresponds to the 
long dashed lines. The short dashed line is obtained when the extra couplings that 
appear in the Euler equation are dropped. They reproduce the large k behavior. 
The dotted line correspond to the case when only the intrinsic coupling of the 
cp field is preserved. It gives the dominant contribution when A is large (right 
panels) but a negligible one when A is small (left panels). 



where we have introduced 



(k2.k) a^m\(P) 



and 



<SEui.(ki, ka) = -^S{ki)r]{k2) (5.30) 

"■1 "'2 



5l„tr.(kl,k2) = ^^^5(fc)77(fc2) 



with 



and 



+ 1.2 S{k^)S{k2)Kk) (5.31) 

"■1 "'2 

The parameters S{k), r]{k) and n{k) have been chosen so that they remain finite when 
k ~ kc- We can now compare the Euler equation derived here with the ones obtained 
in GR. There are two main effects which signal the modification of gravity. The first 
one is the appearance of the e term which was already captured in the 7-model and the 
linear approximation. In the context of scalar-tensor theories, we find that there are 
new contributing terms which are intrinsically due to the non-linearity of the theory. 
Notice that at the second order level, the new terms in the perturbation equations 
amounts to a modification of /3(ki,k2): 

/?cff(ki, ka) = /3(ki, ka) + 5Eui.(ki, ka) + 5intr.(ki, ka) (5.34) 

These terms originating from the new (5m^)^ contributions lead to new effects which 
will be analysed for dilaton models in the following section. 
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5.^. Application to dilaton models 

The dilaton models are good benchmarks to test the influence of modified gravity at 
second order in perturbation theory. To do so, the mass and couplings must be known 
in the cosmological background, it turns out that the mass of the dilaton satisfies [24] 

3^2 + 4»A 

m " 2 (A-2+3(ii^+4)-) ^'-'"^ 

while 

This allows one to deduce the other functions which appear in the perturbation 
equations. Assuming that both A2 and are large, we have 

S(k) ^/^^ ^hen k/kc 00 (5.37) 

They are shown in the various panels of Fig. [5] for A = 10. It is to be noted that S{k) 
is finite in the modified gravity domain. It takes the asymptotic value of about 0.2 for 
VQ = 0- 

The results of the second order modification of gravity are to a large extent 
determined by the asymptotic behaviour of the functions ri{k) and /i(fc). Two cases 
can be distinguished. When A is large the /3off parameter becomes </> independent and 
7eff vanishes. In this case we have 

2 4 - 3^]„ 

9A2 (1 - n„ 



^W"^7TT2 n,o '^2 ^(fc) (5-39) 



^(^^-^ ^^„.(4-3».) ^W- ^'-''^ 

On the other hand, when A is small all parameters tend to be small and follow a 
hierarchical behaviour leaving the terms containing rj as the main coupling effect. In 
this case we have 

^ 3r>^(i - 1 

- - ^ 2(4 -3an) l + a2m2(^)/fc2 ^^"^^^ 

— ^ a.(4-3L)- ^(^)- ^'-''^ 

All the extra couplings actually vanish in the modified gravity regime {k ^ kc) because 
of the a'^m^/k^ factors that enter in the expression of the coupling functions. 

It is also worth noting that the extra coupling functions iSintr. and 5eu1. all vanish 
when ki + k2 = 0. In these cases the modifications of gravity reduce to the ones of 
the linear theory. As a result the effect of the non-linearity of modified gravity are 
only present in an intermediate regime around kc- 

The resulting shape of the bispectra are presented in details in the following (in 
particular in the appendix). The benchmark model we have adopted corresponds to 
A = 10. The value of A2 determines the value of kc and is a free parameter. Most of 
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the results we present are actually independent on the value of k^.- Nevertheless kc 
determines the slope of the primordial density spectrum and so affects the detailed 
results of the bispectra (for non-equilateral configurations). The bispectra are obtained 
as in previous section. For the velocity divergence it corresponds to the reduced 
velocity divergence 9{k)/fk- We can see that although the modified gravity regime 
changes the linear growth rate of the density contrast as soon as fc > fee, the new 
couplings terms play a role only for A; ~ fee for equilateral type configurations (as 
illustrated in Figs. [7|). 
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Figure 8. The density and reduced velocity divergence bispectra in the squeezed 
limit for A = 10. Conventions are the same as for Fig. [3 

In the squeezed limit it is interesting to see that there are significant changes with 
respect to GR in configurations where /cshort S> kc and fciong ^ kc that are due to both 
changes in the linear growth rates and to the extra couplings functions. In Fig. |8]we 
show the case /cghort ~ 10 kc and fciong ~ 0.1 kc- 



6. Conclusion 



We have explored here the consequences of modified gravity models on the evolution 
of mode coupling amplitudes at large scale. In this paper we have concentrated our 
efforts on the computation of the second order expression of the cosmic fluids. For 
Gaussian initial conditions, the induced kernel is expected to determine the amplitude 
and shape of the observed bispectra at large enough scale. 

We found that modifications of gravity can change the amplitude of the 
coupling parameter at different levels. From a pure phenomenological point of view, 
incorporating a change of gravity into a change of the amplitude of the Euler equation 
source term leads to a mild effects in the amplitude of the bispectra. We explicitly 
computed this effect. The results are encoded in Eqs. p.36ll3.37|) in the context of 
what we called the 7-model. 

Realistic models however lead to a richer phenomenology. The change of gravity 
is in general both time and scale dependent. This is the case in particular for the 
chameleon and dilaton models in which the effective strength of gravity is modified 
through extra scalar degrees of freedom in which a scale dependent critical wavemode 
appears which controls the effective gravity strength. It had been stressed in previous 
papers that this scale could be of cosmological relevance. We see here that time 
and scale dependent effects can play a significant role in the amplitude of the mode 
couplings for scales below the critical scales. 

Furthermore the existence of extra interaction fields lead to a more profound 
change of the coupling structure. In particular there exist couplings of arbitrary order. 
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not only quadratic as for standard gravity dynamics at sub-horizon scales (see |37] for 
details.) This is a situation comparable to that encountered in the DGP models and 
described in [3S] . Qualitatively, changing the law of gravity cannot therefore be 
captured in a simple change of the linear theory for scales that are comparable with the 
critical scale. It also implies that the resummation schemes developed in the context 
of standard gravity (see [33l [SU |40l IHl |42] ) cannot be directly applied. 
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Figure 9. The density normalised bispectrum for k2 = 2ki and as a function of 
their relative angle. This is to be compared with the result presented in [39| . 

Quantitatively we find that the coupling functions F2 and G2 that describe the 
2nd order expression of the cosmic fields can be changed at a few percent level 
compared to the standard GR case. The type of changes is however quite model 
dependent and results in the interplay between transient effects and direct coupling 
effects. In Fig. [S]we present the relative change in the density bispectrum for a specific 
set of configurations. These results are be directly compared to the one found in the 
context of the DGP model. We can see that their amplitude is comparable although 
they differ in detailed scale and shape dependences. 

Whether such effects are actually observable in current or future surveys is largely 
an open question. Measuring bispectra, with assessed error bars, proves to be difficult 
and very little has been done or attempted so far. One example is to be found in [43J 
where the bispectrum is measured in the PSCz survey. It has not yet been measured 
in the SDSS as such. Only the angular space three-point function has been measured 
in [44' from which one can infer the error amplitude of such measurements. They are 
in the percent level range. This gives us hope that large-scale surveys in preparation 
can potentially shed lights on these models. 
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Appendix A. The bispectra, full dependence 



In this appendix we simply present the full wavemode dependence of the density and 
reduced velocity divergence bispectra when compared to the general relativity case. 
The modified gravity model we adopted corresponds to A = 10. We assume here that 
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the pivot term for k is and we assume that at that scale the power spectrum is 
close to a power law of index n = —1.5. We do not expect though that the results are 
strongly effected by such an hypothesis. 
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Figure Al. Amplitude of the reduced bispectrum Qs{ki, k2, k'^) for the density 
field for the modified gravity model (A = 10) divided by the expected result for 
standard gravity as a function of k\/k2 and k-^/k2 for respectively k2 = .Ifcci 
k2 = kc and k2 = lOfcc- We assume here that P{k) ~ k~^'^. 

The results are presented in Figs. lAll and |A2I which show the ratio between the 
modified gravity regime and standard GR. Deviation from GR are observed to be at 
the percent level for the density field; up to a few percent level for the reduced velocity 
field. Note that observed density fields in redshift space correspond to a mixing of 
these 2 contributions (sec for instance [57]). 
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Figure A2. Same as previous figure for the reduced velocity divergence. 
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